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Abstract 

As a consequence of QCD factorization theorems, a wide variety 
of inclusive and exclusive cross sections may be formulated in terms 
of a universal colour dipole cross section at small x. It is well known 
that for small transverse size dipoles this cross section is related to 
the leading-log gluon density. Using the measured pion-proton cross 
section as a guide, we suggest a reasonable extrapolation of the dipole 
cross section to the large transverse size region. We point out that the 
observed magnitude and small x rise of the gluon density from con- 
ventional fits implies that the DGLAP approximation has a restricted 
region of applicability. We found that 'higher twist' or unitarity cor- 
rections are required in, or close to, the HERA kinematic region, even 
for small 'perturbative' dipoles for scattering at central impact param- 
eters. This means that the usual perturbative leading twist descrip- 
tion, for moderate virtualities, 1 < Q 2 < 10 GeV 2 , has rather large 
'higher twist' corrections at small x. In addition, for these virtualities, 
we also find sizeable contributions from large non-perturbative dipoles 
(b > 0.4 fm) to i<2, and also to F^. This also leads to deviations 
from the standard leading twist DGLAP results, at small x and mod- 
erate Q 2 . Our model also describes the low Q 2 data very well without 
any further tuning. We generalize the Gribov unitarity limit for the 
structure functions of a hadron target to account for the blackening of 
the interaction at central impact parameters and to include scattering 
at peripheral impact parameters which dominate at extremely large 
energies. 



1 Introduction 



In the proton rest frame, at small enough x = Q 2 /W 2 and Q 2 3> m 2 , Deep 
Inelastic Scattering of a virtual photon from a proton may be viewed as 
being factorized into a three stage process: the formation of a state which in 
general is build of quark, antiquark and gluons from the virtual photon, the 
scattering of this state off the static proton and the subsequent formation 
of the hadronic final state. In QCD, the different quark-antiquark-gluon 
configurations in the photon clearly have different interaction strengths with 
a target. 

For states of very small transverse spatial size, b 2 , the dominant scatter- 
ing state is a quark-antiquark colour dipole state (in this simple case b is the 
transverse diameter of the dipole). This small dipole has a small scattering 
probability, which has been calculated in perturbative QCD, and at small x 
is related to the gluonic colour field associated with the bound state. 

The fact that such a dipole cross section appears in a wide variety of hard 
small x inclusive, exclusive (e.g. heavy vector meson production, DVCS, etc) 
and diffractive processes is a consequence of the feasibility of separating 
scales in QCD [|IJ, and may be formulated in terms of a universal dipole 
cross section at small x. So far, general, relatively ad hoc, ansatze have 
been ascribed to this quantity [||, |2| and the phenomenological parameters 
specified by successful fits to structure function data. 

States of larger transverse size will in general have much larger cross 
sections and will contain many constituents. With the increase of the size 
of the quark-gluon configurations the number of degrees of freedom in the 
photon wavefunction is also increasing (becoming very large in the non- 
perturbative QCD regime as a consequence of spontaneously broken chiral 
symmetry in QCD). Nevertheless, the transverse size of the scattering state 
seems to be an appropriate parameter for the smooth matching between 
cross sections in the soft and hard regimes. For convenience we will continue 
to refer to b 2 as the dipole size, and the cross section as the dipole cross 
section, but these terms should be understood to refer to the transverse size 
and cross section for more general scattering systems for the case of large 
systems. 

The aim of this paper is to exploit the QCD relationship between the 
small dipole cross section (DCS) and the gluon density in the proton to 
build a realistic ansatz, monotonically increasing in b, for the DCS of all 
transverse sizes. In the large b region, we match onto the measured pion- 
proton cross section (at b n = 0.65 fm), to which we ascribe a gentle rise with 
energy. We suggest a smooth interpolation in b for the DCS from small b 
(where perturbative QCD is valid) up to b n and a smooth extrapolation for 

*In practice, the kinematical effect of skewedness of the amplitude in exclusive processes 
leads to a partial, but controlled, breakdown of this universality. 
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even larger transverse sizes. 

Using our ansatz, we then analyze the small x structure functions, cal- 
culated in 6-space which we denote and F%. In this picture, they are 
given by convolutions in b 2 of the square of the light-cone wavefunction of 
the virtual photons, of the appropriate polarization state, with the DCS. We 
observe that for large Q 2 , reproduces approximately the same result as 
the standard leading-log perturbative QCD formula. This agreement pro- 
vides a justification for the relationship between four- momentum scales for 
the gluon density and transverse dipole sizes which we assume in our ansatz 
(Q 2 = W/b 2 , cf. |2], H|, ||). We then produce values for F2 using our ansatz 
and, without any fitting, find reasonable agreement with the HERA data at 
small x, even for the region of low photon virtuality below the input scale 
for QCD evolution, Q 2 < Qq, where our ansatz may be expected to do less 
well. 

It is straightforward to calculate how much of the non-perturbative re- 
gion at large b contributes to the structure functions. Since F<i is mainly 
governed by the transversely polarized photon, the spin structure of the "ij-qq 
vertex leads to a considerably broader integral in b than in longitudinal case. 
We illustrate this well known effect graphically (see figs. (1C, 11, 12])). We find 
that for relatively high Q 2 = 4 — 10 GeV 2 a surprisingly big contribution 
to the integral (as much as 50%) is coming from this 'dangerous' large b 
region. While much of this non-perturbative piece is attributable to the 
input distributions in F2 , the fact that it is also present in Fl could indicate 
a sizeable 'higher twist' contribution for these virtualities coming from the 
non-perturbative region. This presents a severe challenge to the use of low, 
and ^-independent, input scales in the conventional parton density fits (e.g. 
the recent MRST analysis |?j uses input scale Qq = 1.0 GeV). 

We note that the small x rise of the structure function F2 observed at 
HERA, which may be translated into a large and steeply-rising gluon density 
at small x, quickly (after only a few gluon radiations in the ladder) leads to 
a contradiction with unitarity because the DCS for small dipoles becomes 
of the same size as the pion-proton cross section (which grows much more 
slowly with increasing energy). To avoid this problem, it is necessary to 
tame the small x growth of the perturbative DCS. We suggest a way of doing 
this which modifies our ansatz for the DCS at small x. For small enough 
x, this taming appears to be required within the weak coupling limit and 
hence may be related (albeit indirectly and within a restricted range of x) to 
the four-gluon to two-gluon calculations of Bartels and collaborators (for a 
discussion of perturbative higher twist effects in QCD, recent developments 
and references see ||). 

The unitarity limit for the cross section of a spatially-small colourless 
wave packet with a hadron target (see section(8) of Q), and especially the 
theoretical analysis of the amount of diffraction in the gluon channel pre- 
dicted by QCD [10 1, show that the restriction to the leading power of l/Q 2 
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should break down for small enough x, possibly within the kinematics of 
HERA. Moreover, it follows from the application of Abramovsky-Gribov- 
Kancheli [11] cutting rules that accounting for the next power in 1/Q 2 (as 



in iJ) can lead at most to a 25% reduction of the leading twist result without 
introducing negative cross sections j|. Thus, it appears that the decompo- 
sition into leading and higher twists becomes ineffective in the kinematical 
region which can be achieved at the next generation of proton accelerators 
(LHC) or maybe even at the edge of the kinematics of HERA. 

The use of the optical theorem for the scattering of small size wave 
packets off a hadron target makes it possible to deduce a limit (which is an 
analog of the Froissart limit for hadron-hadron scattering) for the amplitude 
in DIS and to calculate the boundary for the applicability of perturbative 
QCD in small x region Q. It was found that for x ~ 10 -4 , the boundary is 
Q 2 ^ 10 GeV 2 . This estimate suggests a significant contribution from higher 
twist effects in the kinematics of HERA for x < 10~ 3 . A more general aim 
of this paper is to visualize this problem and to evaluate structure functions 
of DIS at very small x. We show that many features of the very small x 
behaviour of structure functions can be understood in terms of the geometry 
of the spacetime evolution of high energy QCD processes. 



In the black limit approximation, Gribov [12] deduced the following for- 



mulae for the unitarity limit for structure functions of DIS[|: 

2irr 2 N Q 2 f Ss M 2 dM 2 



Here is the radius of the nucleon and p is the normalized spectrum of pro- 
duced hadronic masses : p(M 2 ) = a(e + e~ — ► hadrons)/c(e + e _ — ► p + p~). 
The upper limit on the M 2 -integral, which imposes the experimentally- 

(M 2 +Q 2 ) 2 m 2 r 2 

observed sharp diffractive peak: —t m i n B£, 2 ^ 1> leads 

to a generic logarithmic energy dependence {Bjy is the usual diffractive 
slope parameter, mjy is the nucleon mass). Strictly speaking this formula 
is valid for M 2 <C s, or <5 < 1, satisfying the condition that the interaction 
for a hadronic system of mass M is close to the unitarity limit. A similar 



formulae has been obtained for Fl [12] and for the gluon distribution [10] 



It is reasonable to ask if, and if so at which x, the black limit will 
begin to be approached in Deep Inelastic Scattering of a virtual photon 
off a proton. In other words, at which point does the cross section for 
the scattering of a small colour dipole with the proton target move away 
from being transparent (due to colour screening) and start to blacken to 
its geometrical limit (a to t = cr d + a ind = 2a inei = 2n(r N + b) 2 fs 40 - 

' Gribov considered scattering off a heavy nuclei for which the black body limit appears 
more natural than for a nucleon. However, provided one assumes the black limit, Gribov's 
arguments, and hence the formulae, will hold. For a recent discussion of the black body 
limit in QCD for DIS off heavy nuclei see flOl n3] and references. 
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50 mb) ? This question is especially acute for Fl and Si^/dlnQ 2 where 
the interaction of spatially small configurations in the wavefunction of the 
photon dominate. We aim to address this question in a phenomenological 
fashion in this paper and to generalize Gribov's unitarity limit to QCD by 
accounting for QCD phenomena which are neglected within the black body 
limit (see section(6)). 

This subject has a rich and long history. For the high energy scattering 
of a hadron from a nuclear target many configurations in the wave function 
of the fast hadron contribute and it is convenient to characterize the interac- 
tion by a distribution, P(cr), of scattering probabilities, a, of its constituent 
states instead of by the average value of a (this useful realization pre-dates 
QCD, see [14]). The qualitative idea of two-gluon exchange as the mediator 
was suggested by both Low and Nussinov [HJ]. Low [15] also observed 
that the dipole cross section should be proportional to the transverse area 
of the object. Miettenen and Pumplin [17] later suggested that scattering 
eigenstates should be identified with partonic configurations in the scatter- 
ing systems, implying that the scattering cross sections of particular states 
should be related to parton densities in the opposing hadrons. In the modern 
context, for DIS at small x and for sufficiently small 6 2 , a oc b 2 a s xg(x,b 2 ) 
and P(c) follows unambiguously from the QCD factorization theorem for 
hard exclusive processes Q. For large b 2 , our approach is in many respects 
similar to the aligned jet model of [18] or QCD-improved aligned jet model 
of | IS] where the cross section of small x processes in the non-perturbative 
regime is expressed in terms of universal dipole cross section at small x, 
which is matched to the soft meson-nucleon cross section (the similarity 
holds despite the difference in the source of qq pairs). 

The paper is organized as follows. In section(2) we discuss the structure 
function Fl in 6-space, introducing a toy model for the DCS for illustrative 
purposes. Section(3) sets out our realistic ansatz for the DCS in detail, 
section(4) compares and contrasts it to other models and ideas in the lit- 
erature. We make some specific, reasonable choices concerning some of the 
uncertainties involved in specifying the DCS. These choices are necessary 
in order to make quantitative statements. However, we have also analyzed 
the precise form of the DCS in detail numerically and investigated the sen- 
sitivity of our results to different choices. This analysis will be presented 
in a separate paper 2C]. However, at this point we merely state that the 
qualitative statements that we make about unitarity and the influence of 
large dipoles in Fi at small x are robust with respect to changes in the de- 
tails. We discuss this in more detail in section(6), where we also consider 
the kinematic region in which the Gribov's black limit may be reached for 
scattering at central impact parameters in DIS for some configurations in 
the photon wave function. We point out that certain diffractive processes, 
for example exclusive photoproduction of J ftp, may act as useful precursors 
to the onset of this new QCD regime. We conclude in section(7). 
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2 Basic Formulae and a toy ansatz for the dipole 
cross section 



In this section we examine the cross section <tl(x, Q 2 ) in 6-space using a very 
simple toy model for the extrapolation of the DCS to large b. Our aim is 
to familiarize the reader with the 6-space formulation of structure functions. 
For clarity of presentation we employ a very simple and unphysical ansatz 
for the DCS at large b. We will improve on this toy ansatz in the next 
section. 

It is convenient to use the impact factor (6-space) representation first 
introduced by Cheng and Wu in considering high energy processes in QED. 
In 6-space the longitudinal structure function, Fl(x,Q 2 ), may be written 
[21, 22] in terms of the DCS convoluted with the light-cone wavefunction of 
the virtual photon squared: 

a L (x,Q 2 ) = 2 I dz d 2 ba(b 2 )\^ L (z,b)\ 2 , (2) 



o 



where 



n.f 

\^(z, b)\ 2 = ^«e. m . £ e 2 q Q 2 z 2 (l - z) 2 K 2 (eb) , (3) 

9=1 

in which e 2 = Q 2 (z(l — z)) + m 2 and for now we set the light quark mass, 
m q to zero. 



For small dipoles, the DCS is governed by perturbative QCD [23 
(for an explicit derivation see |2~ 



7T 2 - 

v pq cd(b 2 ,x) = y b 2 a s (Q 2 ) xg(x, Q 2 ) . (4) 

We employ a phenomenological scaling ansatz Q 2 = A/6 2 to relate transverse 
sizes to four-momentum scales (it is possible to prove that this ansatz is a 
property of the Fourier transform in the LO and NLO approximations but 
not beyond). We also implicitly assume that the DCS is independent of 
light-cone momentum sharing variable z. This is a good approximation for 
the longitudinal case because the average z ~ 1/2 dominate in the integral 
and due to the z — > 1 — z symmetry of the wave function. For the transverse 
case the end points give a larger contribution and hence this assumption is 
less justified. 

The relationship of eqs.(H|||]) holds to leading-log accuracy in Q 2 (so, 
for consistency one is forced to use only LO partons and a s at one loop) 
and involves taking the imaginary part of the usual box and crossed box 
graphs. As such this form corresponds only to the dominant inelastic piece 
of the DCS. We immediately see a practical problem using eq.M) in the 
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b— integral of eq.(^). There are always regions in the integral, at large b, 
where the gluon density is not defined and we need to decide what to do. In 
particular, for fixed A, the gluon density is not defined for b 2 > 6q = A/Qq. 
In the usual treatment this contribution is absorbed into the initial condition 
of the evolution equations. 

To get started we fix A = 10 and simply freeze a s (Q 2 ) xg(x,Q 2 ) at its 
value at Qq for b 2 > 6q in eq.(||) (we refer to this as ansatz 1). This means 
that the DCS retains the canonical 6 2 -behaviour at large b in eq.(||), and its 
derivative is discontinuous at b = 6qo- Fig-(H) shows a plot of the resultant 
DCS function of b for several values of x. 




Figure 1: Dipole cross section in mb for fixed A = 10, with the toy model 
ansatz at large b. For an input scale of Qq = 1.6 GeV, 6qo = 0.39 fm marks 
the boundary of the perturbative region. 



We may now examine the dominant regions of the b- integral in eq.(H): 

poo 

a L (x,Q 2 )= / dbI L (b,x,Q 2 ), (5) 
Jo 

where, 

I L (b) = 2vr6«T(6 2 )/ 7 , L , (6) 

in which I 7) £ is the integral of \ipL(z,b)\ 2 over z. Fig.@ shows I 7i l as a 
function of b for three light flavours at fixed values of Q 2 , it diverges at small 
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values of b due to the logarithmic divergence of the Kq Bessel function at 
small values of its argument. 
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Figure 2: Longitudinal photon wavefunction squared, integrated over z (in 
units of fm" 2 ), for Q 2 = 4, 10, 40 GeV 2 . 



Figs.(^Q) show the integrand, for two characteristic values of Q 2 = 

4,40 GeV 2 using CTEQ4L gluons ||§. Note m fig-© one can clearly see 
that the region above 6qo ~ * 0.197 ~ 0.4 fm (where there is a kink due 
to the freezing of ansatz 1) contributes significantly to the whole integral. 
In contrast, fig.(B) shows that for Q 2 = 40 GeV 2 (typical effective scale for 
T photoproduction ||) this region is completely irrelevant. This is due to 
the fact that the photon piece of the integrand multiplying the DCS, 
strongly weights the integrand to progressively smaller b as Q 2 increases (see 

fig.®)- 

As an example, let us focus on the case when Q 2 = 40 GeV 2 and x = 
10 -3 . The b- integrand exhibits a strong, slightly asymmetric peak around 
small b = frp ea k ~ 0.08 fm, which is slightly skewed to larger b. The 



relationship Q 2 = A/6 2 , with A = 10, implies Qp ea jj. ~ 60 GeV 2 , and that 
the typical b = 6^yp 0.1 fm corresponds to Qjjyp = 40 GeV 2 . Clearly, for 
b < &typ ^ ne effective scale will be larger than Q 2 . The fact that there is 
very little contribution from the large b-region, for large Q 2 , illustrates QCD 
factorization in 6-space: the sharply peaked photon piece of the integrand 



8 



3 



0.012 

0.01 - 
0.008 - 
0.006 - 
0.004 - 
0.002 




"i 1 1 1 1 1 1 1 r 



= 4, x = 


io- 2 


x = 


10" 3 


x = 


10" 4 


x = 


10" 5 



J I i i \-^-—>-— ~-y- - - 1 - - - - - 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

b (fm) 



Figure 3: Integrand of ctl, in units of fm, for fixed A = 10, Q 2 = 4 GeV 2 
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Figure 4: Integrand of cr^, in units of fm, for fixed A = 10, Q 2 = 40 GeV 2 
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ensures that only small dipoles contribute significantly in the integral. 

We will refer to the structure function Fl(x,Q 2 ) calculated in 6-space 
as Fi. It is related to the defined cross section in the following simple way: 



F L (x,Q 2 



F b L (x,Q 2 ,X) 



Q 2 



4vr 2 



-a L (x,Q 2 ,X) 



(7) 



where we have chosen to write the dependence on A explicitly. Why have 
we chosen A = 10 ? Lambda should reflect the typical size of contributing 
dipoles. It may be calculated by defining an average b in the integrand for 
Fl (in [||, H ||] a median average of the integral was used). Whatever the 
precise procedure used for defining the average, a value of A 10 comes out 
for large enough photon virtuality Q 2 . Roughly speaking A =<b> 2 Q 2 , so 
that when b =<b> the gluon and a s are sampled at Q 2 in eq.(||) as in the 
usual leading-log in Q 2 perturbative QCD formula for Fl (see eq.(S) later). 
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Figure 5: The function f(Q 2 ) = a s xg at fixed x for various parton sets 



In fact Fi has a rather weak dependence on A for large Q 2 , this fact 
reflects the renormalization group invariance of QCD. We illustrate this in 
fig.(H) by plotting a s xg as a function of its argument at fixed x = 10~ 3 for 
several LO and NLO parton sets [[?], pq ]. A similar behaviour is observed at 
other small values of x < 10~ 2 . For large Q 2 it is rather a weak function, 
and our A-ansatz translates this into a weak dependence on b. This in 
turn implies that the canonical a pqc d oc b 2 a s (\ 2 /b 2 ) xg(x,b 2 ) oc b 2 holds 
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approximately for small enough dipoles. For x < 10~ 3 , a s xg has positive 
scaling violations implying an effective behaviour slightly softer than b 2 from 
QCD at sufficiently small x (the effective Q 2 -powei, 7(2;, Q 2 ) < 0.1 in fig.(g), 
for large Q 2 > 20 GeV 2 ). The peak in b evident in fig.(|j) implies that a s xg 
is sampled dominantly only in a small range of Q 2 in fig.(||). For small 
dipoles a reasonable change in A corresponds to a shift of this dominant 
region within the fairly flat part of the curve. Hence, as we have checked 
explicitly, for large Q 2 Fl is insensitive to the precise choice of A. 

However, for smaller scales a s xg has a much stronger dependence on Q 2 
which tames this linear dependence on b 2 into something much softer. This 
fact is apparent in the shape of the curves in fig.(|l|), for dipole sizes corre- 
sponding to the region of Qq < Q 2 < 10 GeV 2 , a pqc d deviates considerably 
from the b 2 behaviour apparent for small (b <C &qo) and large (b > 6qo) 
dipoles in our toy ansatz. In [20], we investigate this question of interrela- 
tion of these scales in more detail. Finally, we note in passing that a s xg, 
plotted in fig.([|), has a comparable numerical value and shape for both LO 
and NLO gluons. Since the dipole cross section is proportional to this quan- 
tity for small b, all statements that we make about the size of the DCS using 
LO gluon densities also hold at NLO. 



3 A realistic ansatz for the DCS 

At small x, because the gluon density dominates over the quark density, to a 
good approximation the LO perturbative QCD formula for Fl(x, Q 2 ) which 
we denote generically as F^ is: 



F L (x,Q 2 ) = 4asi f )TR Zf = {e 2 q f dx'x'g(x',Q 2 )^{l-^). (8) 

ZSK J x XX 

This conventional expression for Fl(x,Q 2 ) involves an integral over the 
gluon momentum fraction x , of the proton's momentum p, which feeds 
into the quark box. In fig.@ we plot the integrand of the above formula 
versus x'/x, at x = 10~ 3 , for various Q 2 values. The gluon is sampled at 
a range of values of x > x with the integrand, /|(V/x), peaked around 
x 'peak ~ 1-3 x i an d skewed to x' > x pea k due to the factor multiplying the 
rising gluon density. We define an average, <x>, to be that x 1 up to which 
one must integrate to obtain half of the full integral: it turns out this is 
always around x =<x>~ 1.75 x for a wide range of external x, Q 2 . 

Momentum conservation for the fusion of a gluon with a photon, of 
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momentum q, to produce the quark-antiquark pair, of mass M 2 g , gives 



,2 i ™2 



kf + m 2 



(x'p + qf = M 2 qq 

z I 1 — z ) 

k 2 I 



z(l-z) 



The inequality in the second line is satisfied for a non-zero quark mass, 
m q , when z = 0.5. In fact the approximation l/(z(l — z)) ~ 4 holds over 
a reasonable wide range of z values. So, taking finite quark masses into 
account implies a minimum value for x of x min = x (1 + 4m 2 /(J 2 )- To 
account for the important role at large b of confinement and spontaneously 
broken chiral symmetry we choose constituent quark masses. With this in 
mind, for our new ansatz we choose to sample the gluon density in eq.(||) at 



x ' = x ' min (l + 0.75^^). (10) 

This choice guarantees that for the average b 2 , x' =<x>, but allows x (b 2 ) to 
vary according to the inverse of the transverse size of the dipoles. Kinemat- 
ically, a large mass dipole requires a gluon carrying a greater than average 
momentum fraction to produce it and the ansatz of eq.(|To|) is designed to 
reflect this. In contrast, for very large, small mass dipoles M 2 q <C Q 2 , and 
the formula gives x' ~ x' min , which is approximately x only for light quarks. 
The average b is <b> 2 = \/{Q 2 + 4m 2 ,), in agreement with the definition 
of A =< b> 2 Q 2 in the DIS region, Q 2 3> 4m 2 ,, but is infra-red safe in the 
photoproduction region Q 2 <Am 2 q . 

We also have a constraint on the DCS at large b from the experimental 
determination of the pion-proton cross section, = 23.78 mb ptSfl . The 
DCS should be close to this at transverse separations which correspond to 
the diameter of the pion (d^ ~ 0.65 fm). For rather large x ~ 10 -2 the 
magnitude of the DCS at the interface to the non-perturbative region is 
considerably smaller (a pqc d(x = 10~ 2 ,6q ) 6 mb) than this. 

The intermediate region 6qo < b < b^ is extremely interesting and very 
poorly understood, so some modeling is required. It is in this region that 
strong confinement and effects of spontaneously broken chiral symmetry set 
in to produce the bound state pion. Clearly the dynamics in this region 
will include strong colour fields and the creation of light sea pairs from the 
vacuum. As such, it is no longer reasonable to think of b as corresponding 
to the transverse size of a dipole. It is better to think of it as corresponding 
to the typical transverse radius of the complicated non-perturbative system, 
which in general will contain many constituents. 

The minimum requirement of a interpolating function, &j for the DCS 



12 



2000 



1800 - 

1600 - 

^ 1400 - 
1200 



1000 



■JS 800 
^ 600 



t: r 

Ql = 45 - 
Q 2 = 10 - 
Q 2 = 4 - 




Figure 6: Integrand of F[ versus x jx for various Q 2 using CTEQ4L gluons. 



is that it matches appropriately at b = 6qo an d b = b n : 

°l(x, b 2 ) = [(Tnp{x, bl) - <7 pqc d( & Qo)]^( &2 ) + ^pqcd( 6 Qo) i ( n ) 

with H(bQ ) = and H(b 2 ) = 1. On geometrical grounds, we also choose 
to only consider functions which are monotonically increasing as a function 
of b. A very simple function which satisfies these requirements is 



Hi(b 2 ) 



(b 2 ~ b 2 Q0 ) 
(bl ~ b 2 Q0 ) ' 



(12) 



which has a linear growth in b 2 , even for b ~ b w . To impose a natter 
behaviour in this region and a fairly smooth matching close to b ~ 6qo we 
choose the following exponential matching 



H( 



e-1 



-[l-eM-Hi(b 2 ))}, 



(13) 



which retains the linear growth in b 2 close to 6qo- 

The pion-proton cross sections is observed to rise slowly as the energy 
increases. In order to take this into account we impose a slow growth with 



increasing energy, consistent with a Donnachie-Landshoff soft Pomeron 1 27 ] 
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in our boundary condition at b = b 1T : 

O.ON 

X J 

and choose xq = 0.01. This behaviour in x is designed to mimic the observed 
(W 2 /Wg ) 0-08 behaviour in energy. The precise value of 23.78 mb is taken 



(T^,n(K,x) = 23.78 (-2) mb, (14) 



from the Fermilab data 2J] from pion-proton scattering and corresponds to 



Wl = 400 GeV 2 . 

As x decreases the magnitude of the DCS for small b 2 increases much 
more rapidly with energy than this soft piece, as a result of the steeply ris- 
ing gluon density. For leading-log gluons at small enough x, if unchecked, 
it can even become greater than ^^(b^, x) for perturbative b < 6qo- This 
is clearly nonsensical, some taming of this rapid growth must occur before 
this can happen. Since the gluon density form of eq.@ really represents 
only the inelastic part of the dipole cross section, as it becomes comparable 
to cx^jv we need to include the elastic part of the DCS too, which hitherto 
was implicitly assumed to be negligible. In the limit of very large energy 
the Froissart bound indicates that the elastic piece should not exceed the 
inelastic piece. On this basis we really should worry about the applicabil- 
ity of our perturbative QCD formula when the DCS is about 50% of the 
pion-proton cross section. The absolute upper bound for the inelastic small 
dipole-nucleon interaction is 8irB2 g /(l + rf) where B2 g ~ 4 — 5 GeV~ 2 is 
the slope of the i-dependence of the two-gluon form factor of the nucleon, as 
measured in hard exclusive diffractive processes at HERA, and ij is the ratio 



of real and imaginary parts of the scattering amplitude fl|, |10[. This bound 
is slightly weaker than the transition point which we assume, but pushing 
all the way to the absolute limit appears unrealistic. 

In our computer code we test a pqcc i (x,b 2 ) to see if the equality is reached 
in the perturbative region b = b cr u < 6qo> where b cr u is defined implicitly 
by 



a(x,b 2 crU ) =^a s (Q 2 crU )x'g(x',Q 2 crit ) = (15) 



with Q 2 rit = ^/b 2 rit and x' is given by eq.(lC). If so, we use a new interpo- 
lation in the region b cr n < b <b- K : 



°i( b2 > x ) = [ h2 \ a 2 ) CT o- (16) 



Matching at b = b n sets the value of 



Mx) = ^AK,x)(^^-) n . (17) 
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The two remaining parameters, a, n, are chosen to provide a fairly smooth 
matching at b = b cr u. To achieve this we perform a three parameter fit of 
exactly the same form as eq.(|l6|) for a given x in the region just below b cr i t , 
using MINUIT |^8| . We then take the effective power, rim, from this fit, so 
that the interpolating ansatz of eq.(|l6|) has approximately the correct power 
in b 2 at the boundary. 

The last remaining free parameter, the scale a, is then specified by the 
matching at b = b cr it- 

2 6 c 2 „ t (l-(0.5)"^) 
a - —2 . 

(^(0.5)"/** - 1) 

This ensures a fairly smooth behaviour in b 2 which takes into account the 
effective behaviour of a s (Q 2 )xg(x, Q 2 ) close to b cr u as discussed earlier (see 

fig-(D)- 

For very large dipole sizes, b > b n , we simply impose a universal residual 
slow growth, linked to the value at b = b n of the form 

aK6 2 >6^) = ^,x) (62 1 + 5 | /2) . (19) 

Numerically this very large b region is totally irrelevant for the calculation 
of DIS structure functions since it is killed by the exponential fall-off of the 
photon part of the integrand due to the Kq, K\ Bessel functions (in practice, 
for moderate Q 2 we integrate up to b = 1.0 fm, for smaller Q 2 < 3.0 GeV 2 
we extend this out to 2.0 fm). 

Let us briefly summarize our realistic ansatz. Assuming the universal 
scaling relation A = b 2 Q 2 restricts the region of applicability of any per- 
turbative QCD dipole formula to transverse sizes smaller than 6qo, which 
corresponds to the input scale of parton densities in 6-space. For small, 
but not too small x, we use the perturbative QCD formula of eq.(Q) with 
the gluon sampled at x (see eq. (JlOj) ) in the region < b < 6qo- Between 
6qo ~ 0.4 fm and b w = 0.65 fm we use the interpolating formula of eq. (|Tl|) , 
employing the exponential matching of eq.(13). For smaller x, we recognize 



that dpqcd gets too large within the perturbative region b < bQo at some 
point, b cr u, defined by eq. (|l5|) . At this point, we switch from the standard 



form and use the interpolating form of eqs.(16,17. 15), which tames the rapid 
growth and interpolates in the region b cr u < b < b^. 

In both cases, our ansatz matches onto the pion-proton cross section of 
eq.(|J), at b = b^ = 0.65 fm, which is allowed a slow Donnachie-Landshoff 
type energy growth. For b > b^ we use the slow increase given in eq.(^). In 
fig-(@) we show our new ansatz for the DCS as a function of x. Note, we use 
the ansatz for both heavy and light flavours. The NA38 collaboration 2£] 
recently suggested N m 24 ± 5 mb, on the basis of an observed deficit in 
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the number of if) decays to dimuons in Sulphur-Uranium collisions relative to 
well established trends in proton-nucleus collisions. The fact that the large 
cc(2S) bound state has a large interaction cross section with the nucleon 
in these nuclear collisions, as predicted in |3(|, is some justification for our 
flavour blind choice for the DCS in the non-perturbative region (6 ~ 0.6 fm). 




b (fm) 



Figure 7: Dipole cross section in mb for fixed A = 10, with the realistic 
ansatz at large b. For small enough x unitarity corrections are included. 



4 Comparison with other models for the dipole 
cross section 

The unitarity correction at small x, discussed above, is clearly beyond the 
usual DGLAP leading twist analysis and is similar in spirit to the saturating 
ansatz of Wusthoff and Golec-Biernat @ : 

a(x, b 2 ) = <7 (1 - exp[-b 2 Q 2 /4(x/x ) x ]). (20) 

A three parameter fit to the HERA data on DIS with x < 0.01, excluding 
charm and assuming Qo = 1-0 GeV, produced the following values uq = 
23.03 mb, x = 0.0003, A = 0.288 and a reasonable x 2 ■ We are encouraged to 
note that the "saturation" cross section coming from this fit is considerably 
below the lowest value of the black limit {o^t}^^ = 2nr 2 ^ ~ 40 mb). 
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Figure 8: Dipole cross section in mb for Wusthoff-Golec-Biernat saturation 
model g. 



The resultant DCS from eq.(|20|) is plotted in fig.(||). A comparison of 
this figure with our model for a in fig.(^) is shown in fig.@ which shows 
the ratio of our a divided by the Wusthoff Golec-Biernat form, a WGB . It 
reveals a considerably different 6-shape, normalization and energy depen- 
dence. Focusing on exclusive processes, which are particularly sensitive to 
the small b region in which the models differ most, will help to distinguish 
between them. Having emphasized the contrasts, it is worth pointing out 
that the two models share some gross features: an approximate b 2 behaviour 
and strong rise with x at small b, tamed to something much softer at large 
b. Indeed, the critical point at which we apply our unitarity corrections, 
bcriti clearly shifts to smaller b as x decreases, result of the rising gluon 
density. This is similar to the critical line of ||]. 

The form of eq.(|20|) was clearly chosen with simplicity in mind and is 
indeed impressively economical in its number of parameters. We recognize 
this motivation and so do not wish to criticize it too strongly. However, 
we feel this simple form misses several known crucial features which our 
ansatz includes (leading to the differences manifest in fig.(||)). Firstly, the 
small dipole form is known from perturbative QCD (having specified the 
relationship between b and Q 2 , see eq.@). The fact that the gluon density 
may be taken from the global fits allows a careful study of the deviation from 
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Figure 9: Ratio of our dipole cross section with the Wusthoff-Golec-Biernat 
saturation model ||. 



a simple b 2 behaviour inherent in eq.(^). As a result the difference between 
our a and a WGB for the perturbative region of b < 0.4 fm is large and 
strongly energy dependent (cf. fig.(|9])). Using the fitted gluon density also 
allows the precise behaviour in x to be incorporated more correctly than 
a single global power implied by eq.([20|). Secondly, it is well established 
experimentally that soft dipole cross section increases slowly with energy. 
This may be modeled with a small power as in eq .(|i~4| ) or with a logarithm 
according to theoretical prejudice. However, eq.(|20|) has a flat behaviour 
in energy for large dipoles, leading to a rather large and energy dependent 
difference for large b. Thirdly, as discussed above, the region where unitarity 
has to be included is known from the requirement of smallness of the elastic 
cross section as compared to the total cross section. 

The eikonal form assumed in [0], i.e. eq.(f20|), was inspired by an earlier 
work by Gotsman, Levin and Maor |HJ (see also J32|| ) in which the hard cross 
section & l ™f d oc a s xg |23j explicitly appears in the eikonal |. At the same 
time the eikonal approximation has problems in accounting for the generic 
properties of QCD. In particular, the eikonal approximation assumes con- 

* However, we remind the reader that if crl^d is large enough to require eikonaliza- 
tion then a^ qcd will also be so large that in the black limit one has a tot = J d 2 p(2 — 
2exp(—a i p ° q t cd T(p)/2)) within the eikonal approximation. 
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servation of bare particles in the wave function of the photon despite the 
fact that operators of Lorentz boosts do not commute with the operator 
of the number of bare particles. A related problem is that the eikonal ap- 
proximation neglects the energy lost by an energetic particle in the inelastic 
collisions. Hence, the taming of the increase of parton distribution by this 
method strongly overestimates the energy released in inelastic collisions. 

The recent paper of these authors with Naftali [33], appears to be some- 
what orthogonal to their earlier works [31, 32] in that the problem of c*™^ 
getting too large is not addressed. They discuss different regions in the 
mass of the scattering state, in contrast to dipole sizes (cf. eq.(|l])). Roughly 
speaking, large masses correspond to small dipoles (at least for the lowest 
qq Fock state) and are related to the unintegrated gluon density in this hard 
contribution. The small mass region is modeled using a Regge analysis. As 
we have discussed in the current paper (see also |3^|), they also stressed the 
interplay of short and long distance physics to all processes at small x. This 
correspondence between masses and regions in transverse size is known to 
break down for higher order Fock states due to the possibility of large size 
and large mass aligned jet type configurations (see e.g. f35|| ). 

An important point of our analysis, which has been stressed elsewhere 



including in |4, 31, p2j, is that it is the magnitude and increase of the gluon 



density within the DGLAP approximation as x decreases that leads to con- 



flict with unitarity. However, in contrast to [|2|, |31j, |32|] we predict that 
structure functions continue to increase significantly with energy above the 
unitarity limit as a result of important role of peripheral collisions. This 
is because unitarity only restricts the contribution to structure functions 
related to the collisions at central impact parameters (see section(6)). 

Forshaw, Kerley and Shaw Q propose a general, rather ad hoc, ansatz 
for the DCS, which they stress is modeled as a function of b 2 and W 2 (rather 
than x), with a soft Pomeron and a hard Pomeron piece: 

P 2 (h) 

&{W2 ' h) = ° 1 + P 2 (b) {b2W2)XS + fe2Pfe2(6) e M-vtb)(b 2 W 2 ) x * , (21) 

where P s (b) and Ph(b) are polynomials in b. They successfully fit this general 
form to the data. From the point of view of very large dipoles (of the order 
of a meson size) one may think that the DCS should be a function of W 2 
rather than x = Q 2 /W 2 , since it should just depend on the energy of the 
collision (here there is no hard scale with which to specify an 'x'). However 
the presence of a finite x imposes off-forward, or skewed, kinematics which 
cannot be ignored even in the limit of soft interactions. One can see this 
for example in the Aligned Jet Model where the propagator for a transition 
from negative mass-squared (— Q 2 ) to positive mass-squared of the qq pair 
is present for DIS. Hence, in order to achieve the observed approximate 
Bjorken scaling of F2 it was necessary to model the DCS in eq.(^) as a 
function of b 2 W 2 , rather than W 2 alone. It is interesting to note that 
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with our ansatz for the relationship between small dipole sizes and hard 
four-momentum scales, this reduces to an x -dependence for small dipoles 
(b 2 W 2 ) n — > (XW 2 /Q 2 ) n oc x~ n . For small dipoles, in our approach, there 
is an identifiable set of diagrams which contain the fusion of a gluon with 
the virtual photon, hence the DCS can and must depend on the momentum 
fraction x of the incoming gluon. How, and whether, this x— dependence 
becomes a W 2 dependence for large dipoles, and whether it is possible to 
define the DCS as a unique function over the whole range in 6, are interesting 
and open questions which deserve further consideration. We also note that 
the DCS which results from the fit of eq. (|2l|) , has the rather unattractive 
feature at large (fixed) energies that it is not a monotonically increasing 
function of b, it contains a minimum in the region b ~ 0.6 fm for s = 
10 5 GeV 2 (see fig. (6,7) of ||). By design, our ansatz specifically avoids this, 
for fixed x. 

From a radically different point of view, we have arrived at a broadly 
similar picture to Donnachie and Landshoff's two Pomeron model with 
the soft Pomeron becoming of increasing less importance as the hardness 
of the process increases (commonly called 'higher twist'). The approximate 
scaling observed in a s xg in fig.(||) explains why a single power in energy will 
work fairly well for the hard piece. At the same time there are important 
differences, in our picture for intermediate b the power is different from 
either of the Pomeron powers of [36] and description of it by the sum of two 
powers is only approximate and looks rather artificial from the point of view 
of the dipole picture. Even more importantly, due to unitarity effects and 
the important role of peripheral collisions we expect a change of the power 
at higher energies making it closer to the soft Pomeron case. For example, 
in exclusive production of J/ij) the logic of our ansatz suggests that the 
large power observed at HERA will become tamed to a smaller power at 
even higher energies due to unitarity corrections. In contrast, in the two 
Pomeron picture of [36] the harder Pomeron would completely dominate at 
higher energies. 



5 Testing our ansatz: comparison with structure 
functions 

In section(2), we argued that for very large Q 2 the integrand in ol is strongly 
peaked in the perturbative region and has very little influence from non- 
perturbative effects in the large b region. If this is really the case, and our 
ansatz is reasonable, we should be able to reproduce values for the structure 
function Fl(x, Q 2 ) which are in close agreement with the standard leading- 
log perturbative QCD formula of eq.(||). As Q 2 decreases towards the input 
scale Qq we might expect the two formula to deviate since the large b region 
is implicitly excluded from the leading twist perturbative QCD formula. For 
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consistency we use the same parton set in each and to avoid the complexities 
of treating massive charm, in this theoretical cross check, we run with only 
three light flavours of quarks (we also set the three light quark masses to 
zero). 

Table. (|l]) reveals the excellent agreement of the 6-space formula with 
perturbative QCD at large Q 2 and also displays the deviation of the two 
formula for low Q 2 . Also shown, in the last column is the percentage, Rl(%), 
of the b-integral coming from the non-perturbative region above b > 6qo- 
As expected this decreases with increasing Q 2 . We used CTEQ4L parton 
distributions which have an input scale of Qq = 2.56 GeV 2 . 
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Table 1: A comparison of the 6-space formula, F^, using our ansatz for the 
DCS, with the standard perturbative QCD result, F?, for F^ f 3 (x,Q 2 ) as 
a function of x for Q 2 = 4, 10,45 GeV 2 . We used CTEQ4L gluons. 

Having found reasonable agreement with this theoretical cross check of our 
ansatz for the DCS, we proceed to calculate its predictions for F2(x,Q 2 ), 
which we denote, F%. In order to calculate this we need to know the wave- 
function squared for transverse photons: 

|^T(^6)| 2 = ^2a e . m .^ e 2 [(z 2 + (l-z) 2 ) e 2 K 2 (e5)+m 2 ^ 2 (e5)] , (22) 

q=l 

where e 2 = Q 2 (z(l — z)) + m 2 . From now on we will use m q = 300 MeV, 
for u,d,s and m c = 1.5 GeV, for both the longitudinal and transversely 



21 



polarized photon wavefunctions. The small light quark constituent mass 
only affects the structure function seriously in the case of very small Q 2 < 
Qg. It acts as a regulator for the divergence of the Bessel function in the 
photoproduction limit Q 2 — ► 0. 

We are interested in how IV'tI 2 ; integrated over z, weights the dipole 
cross section in 

fl/2 roc roo 

a T (x,Q 2 ) = 2tt2 / dz bdba{b 2 )\^ )T {z,b)\ 2 = I dbI T (b,x,Q 2 ) , 
Jo Jo Jo 

(23) 



where ir(fr) = 27t66"/ 7i t and 

r l/2 -1 V™ 7 P 2 

J 7>T (6) = 2 jf dz e - m £r « [{z 2 + (1 - ,) 2 ) e 2 *T 2 ( e &) + m 2 q K 2 {eb)] . 

(24) 

This integrand, multiplied by the Jacobian factor 2nb, is shown in fig,(|T0|) 
as a function of b for three light flavours. At small b, K\{a) oc 1/a so I 7j t is 
approximately independent of Q 2 and the residual 1/b 2 , apparent in fig.(^0|) 
cancels the b 2 in eq.(f|). As a result, in the perturbative region the transverse 
cross section, which dominates F2, is particularly sensitive to the effective 
behaviour of a s xg in b. For comparison, fig.([Tl|) shows 2irb In contrast, 
in the longitudinal case Kq(o) oc — ln(a) at small a, leaving almost the full 
power of Q 2 evident in eq.@, and effects the approximate 6 2 -behaviour of 
a very little. 

In order to make a comparison of the relative weight that the trans- 
verse photon gives to the DCS, we plot the ratio of J 7; t and I 7i l in fjg(|i~2|). 
This plot clearly shows that the transverse photon provides support over 
a much broader range in b, i.e. both at smaller and larger b, than the 
longitudinal photon. This is especially true at high Q 2 and leads to a 
very broadly peaked integrand It- The peak corresponds approximately 
to Q 2 and smaller dipoles b < b pea k strictly speaking lie outside of the usual 
leading-log approximation (which result from logarithmic integrations in the 
perturbative QCD ladder in kt up to Q 2 ; the variable b is conjugate to the 
kt in the upper quark loop of the ladder in this formulation) . 

The fact that the transverse photon wavefunction squared is so broad in 
6-space reflects Gribov's paradox (or Bjorken's aligned jet model). At large 
Q 2 the large dipoles are being produced by asymmetric splittings (« < 1). 
While such splittings are unlikely, the large dipoles they produce interact 
with a large hadronic cross section. The detailed link between large dipoles 
and asymmetric splittings, and the possibility of a z-dependence in a, re- 
quires further study. Figs. (1^,1^) compare the resultant integrands for the 



transverse and longitudinal photons for large (Q 2 = 40 GeV 2 ) and moder- 
ate (Q 2 = 4 GeV 2 ) values of the photon virtuality, respectively. For the 
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Figure 10: The weight given to dipole cross section by the transversely 
polarized photon as a function of transverse size. 



latter case, one can clearly seen the unitarity corrections starting to affect 
the integrands at the smallest values of x. The relative contributions of 
perturbative and non-perturbative regions are very clear from these figures. 
To calculate F^ ^ we use 

O 2 

F 2 (x,Q 2 ) = F T + F L = -^ (vt + ctl) (25) 

with ctt defined precisely analogously to ol- We also attempt to include 
threshold effects, albeit in a rather crude way, by imposing that the momen- 
tum fraction of the gluon must be sufficient to generate the charm quark 
pair, using a theta function Q(x' — (Q 2 + M 2 C )/W 2 ). This condition is in- 
cluded in our ansatz of eq.(10) and ensures that x' > x min = x (l + 4m 2 /Q 2 ) 
for any value of b. This procedure leads to approximately the correct ratio 
of charm in F%. 

In fig.(|l5|) we compare our results (solid curves) with the 1994 HERA 
data |57j, for the larger Q 2 > Qq = 2.56 GeV 2 data. Our physical ansatz 
is in reasonable agreement with the data at small x for moderate values of 
Q 2 . This gives us faith in some of the choices we made in specifying our 
ansatz for the dipole cross section. Recall that no fitting or minimization 
procedure has been applied to tune the available parameters to give an 
excellent fit, although of course it would be possible to do so. 
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Figure 11: The weight given to dipole cross section by the longitudinally 
polarized photon as a function of transverse size. 



In addition, in fig. ([[(]) for completeness we show an explicit comparison 
of our model with the lower Q 2 < Qq ZEUS BPC |]39| , |40tl and a selection 
of the ZEUS SVX El| data. We could equally have chosen to compare to 
the low Q 2 HI data from 1995 [42], which is binned differently in Q 2 . It is 
interesting to note that our model does a reasonably good job in this low 
virtuality region too, which is most sensitive to our ansatz for the DCS at 
large b. Again recall that no fitting procedure has been applied. 

More recent HI data, over a wider kinematic range than in [42] was 
recently presented [43|. Unfortunately the data tables were not publicly 
available to include in our plots. 

Note that in this region of small Q 2 we clearly cannot use DGLAP evo- 
lution. At the same time, due to the possibility of separating contributions 
of small transverse distances in the wave function of the virtual photon with 
small virtuality we can estimate with rather small uncertainties the short- 
distance contribution which leads to a fast increase of the cross section with 
W 2 at fixed Q 2 . 
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Figure 12: The ratio of weights, transverse divided by longitudinal, given 
to the dipole cross section for the two different photon polarizations. Three 
light flavours are included. 
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Figure 13: Integrands in 6-space (units mb fm ) for Q = 40 GeV . 
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Figure 14: Integrands in 6-space (units mb fm r ) for Q 2 = 4 GeV 2 . 
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F 2 (x,Q ) Comparison at low Q 
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Figure 16: A comparison of the results with the ZEUS data at low Q 2 , 
without any fitting procedures. 
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6 The onset of the new QCD regime 



The small x dipole formulation of high energy processes is attractive in that 
the contributions from perturbative and non-perturbative regions become 
very clear. We have extrapolated the perturbative QCD formula for small 
dipoles into the non-perturbative regime, using the pion-proton cross section 
as a guide. For very small x it was necessary to tame the growth of the small 
dipoles to avoid conflict with unitarity. 

In our ansatz, small dipoles are governed by the steeply rising gluon 
density. In principle all high energy processes contain perturbative and non- 
perturbative contributions, which rise quickly and slowly with increasing 
energy, respectively. Having specified the dipole cross section, the light- 
cone wavefunctions of the external particles pick out a given region in b 
according to their hardness. A reanalysis of fairly soft exclusive processes, 
(e.g. electroproduction of p-mesons, or DVCS) will be crucial in constraining 
the large dipole region further. 

So, where does this leave the conventional DGLAP analyses of inclusive 
structure functions which always produce qualitative good fits to the data? 
In our opinion, the linear nature of the evolution equations mean that the 
global fits are almost guaranteed to work. When new data comes out it 
is always possible to fine tune the functional form of the input, and the 
many input parameters at the starting scale, to reproduce the slow loga- 
rithmic changes in Q 2 in F2. However, when the resulting quark and gluon 
densities produce apparent contradictions it is surely time to extend the 
conventional picture to include other (higher twist) contributions. This is of 
course very far from straightforward in practice although some interesting 
attempts have been made |g, 44]. A good way to establish experimentally 
the important role of higher twist effects is to measure diffraction in the 
gluon channel and/or diffractive charm and beauty production in DIS (see 
[]To| for discussion). 

As well as discussing the transverse size of the scattering system it is 
also useful, and usual, to discuss the typical impact parameter, p, of a given 
configuration on the target (i.e. the conjugate variable to the transverse 
momentum transfer q t ; t = —q 2 ). Although it is standard to discuss the 
impact parameter representation of scattering process (see e.g. the textbook 
by Eden |^5| ) , we have noticed that there is often a confusion between p and 
transverse size b in recent literature. Some useful explicit formulae in this 
regard are given in the recent paper by Gieseke and Qiao [46]. For particular 
processes the typical impact parameters are usually expressed in terms of 
the slope parameter B which specifies the i-dependence (assumed to be 
exponential). 

In hard exclusive diffractive processes in DIS, within the kinematics of 
HERA, the interaction is typically dominated by scattering at central impact 
parameters: a(p 2 ) tx exp(— p 2 /2B), where B « 4.5 GeV -2 . Thus, very ap- 
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proximately, for these processes p ~ V-S ~ 0.6 fm which is significantly less 



than the electromagnetic quadratic radius of a nucleon r t ^ ans ' « y 2r 2 N /2> 
0.84 fm. In contrast, in hadron-hadron collisions peripheral (large p) in- 
teractions play an important, and perhaps dominant, role. For example, 
an analysis of elastic proton-proton collisions at Fermilab collider energy 
range shows that B ~ 17 GeV~ 2 and therefore peripheral p ~ \/2B ~ 1.2 
fm are essential in the total cross section. On this basis we can conclude 
that the spacetime development initiated by a spatially small wave packet of 
quarks and gluons in DIS is different at the achieved energy ranges from the 
pattern known from proton-proton collisions, where soft Pomeron physics 
dominates. 

From the analysis of partial waves in high energy (^/s > 50 GeV) elastic 
pp collisions we know that non-perturbative interactions at central impact 
parameters are close to the black limit and that peripheral collisions play 
a crucial role. The distinctive feature of DIS, in the kinematics of HERA, 
which we understand from the above analysis of slopes, is that the scattering 
at central impact parameters dominates and peripheral scattering, where the 
interaction is far from black in a wide range of small X j IS cl correction. The 
advantage of considering the black regime in DIS is that it is possible to 
make a complete evaluation of structure functions of a proton, (at least 
at not too large Q 2 ~ 1 GeV 2 , i.e. in the regime where nonperturbative 
QCD physics dominates in the structure functions) because in this limit 
all configurations interact with the same cross section. Since geometrically, 
large transverse size dipoles are more likely to correspond to large impact 
parameters, peripheral collisions should play a much bigger role in ut than in 
ctl for DIS since the large b plays a much more significant role (cf. fig. (12)). 

To demonstrate that our analysis is rather general and almost model 
independent it is useful to rewrite the equations of the dipole model in a 
form which accounts for the blackness of interaction for b 2 > 6q, i.e. to 
evaluate the black limit of cross section: 



*tot(ll, T +P^X) = J dzd 2 b a(z,x,b 2 )\rf' T (z,b,Q 2 )\ 2 

= 27rr 2 N J dzd 2 b\rf' T (z,b,Q 2 )\ 2 

- J dzd 2 b[27rr 2 N -a(z,x,b 2 )]\^ T (z,b,Q)\ 2 9(b 2 -b 2 ) 

+ J dzd 2 b a periphera i(z,x,b 2 )\^ T (z,b,Q 2 )\ 2 . (26) 

The first term in the last expression is relevant for central collisions and as- 
sumes blackness for all dipole sizes^J : a(z, x, b 2 ) = 2^7^ (including inelastic 

''It is unclear whether the hard amplitude actually achieves the black limit or if the 
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and elastic contributions) . Alternatively, this black limit of Gribov [T^j may 
be expressed in terms of the structure functions 

2irr% r Ss dM 2 p{M 2 )M 2 Q 2 

J (M 2 + Q 2 ) 2 
2ttt% f 5s dM 2 Q*p(M 2 ) 

12tt3 J (M 2 + Q 2 ) 2 ■ { '> 

As we mentioned in the introduction (cf. eq.((l|)) the integral over the mass, 
M 2 , of the intermediate state leads to a logarithm in energy, or equivalently 
in 5/x: 

F 2 = 2 ^^HS/ X ). (28) 

This estimate is rather close to that obtained within the BFKL approxi- 
mation ([^] and references therein) which we don't explore in this paper. 
Based on the equations for the black limit suggested in this paper, which 
account for the characteristic suppression of the interaction for sufficiently 
small configurations we may identify S as the critical scale x cr n — M 2 /s 
at which unitarity corrections become necessary for a given s and M 2 . An 
analysis of the black limit in QCD, to be published in a separate paper, 
shows that for fixed Q 2 5 decreases as x increases, leading to a reduction 
of the coefficient of ln(l/x) by a factor of 2 — 4 from one for a fixed 5 (de- 
pending on Q 2 ). The size of this reduction depends on the rate of increase 
of xg(x, Q 2 ) with decreasing x and increasing Q 2 . 

For Q 2 of a few GeV 2 , x ~ 1CP 5 using our estimates of the kinematic 
range in which unitarity effects may become important we find for 5: 10~ 3 < 
S ~ x cr it < 10~ 4 . This allows us to estimate the numerical values of F 2 for 
which black limit corrections will be important. For Q 2 = 1 GeV 2 x ~ 10~ 5 , 
and taking a reasonable constant p ~ 2.5, one obtains .pblack ~ 15 _ 3 
This estimate is only a factor of about three bigger than the current HERA 
measurements. 

The second term in eq.(|2^) corrects for the fact that for small b < 60 
this black limit will not yet have been reached, at the x value concerned. 
This is completely general and does not require a detailed knowledge of 
the wavefunction of the photon in the large b 2 region. A hypothesis on the 
blackening of the interaction permits a calculation of the structure functions 
at very small x which otherwise can not be evaluated within the existing 
methods of QCD. The second term accounts for the details of QCD phe- 
nomena: because of colour screening phenomenon and asymptotic freedom, 
which are built into the QCD expression for small dipoles (cf. eq.(Q)), the 
Gribov unitarity limit can not be achieved for all configurations in the wave 

increase of the perturbative QCD amplitude with energy is tamed earlier. Either case 
leads to the unitarity limit for collisions at central impact parameters. 



Ft = 
F L = 
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function of a photon. At any particular x, there will exist configurations 
whose interaction is far from the black limit. So, to evaluate this term it is 
necessary to make some additional assumptions concerning how the black 
limit is reached in practice. This introduces some model dependence. For 
example, one may choose to use our ansatz for a(z,x,b 2 ) with its specific 
assumption of a smooth interpolation of it to unitarity limit. 



The third term of eq.(26) concerns peripheral collisions which will dom- 
inate for any initial configuration in the photon wavefunction at extremely 
small x, due to Gribov diffusion in the parton ladder. The reason for this 
asymptotic dominance is that for a particular transverse size central colli- 
sions will freeze at their black limit, whereas peripheral collisions can con- 
tinue to grow (albeit slowly) with energy. For example, for a Donnachie- 



Landshoff soft Pomeron parameterization [27] they would continue to grow 

like Cp er ipheral ^ ( s / s black) e - 

In this paper we have not tried to elaborate in detail the behaviour of 

the structure functions in the kinematics where the increase with energy 

of the perturbative QCD amplitude is slowed down (i.e. where the dipole- 

nucleon cross section is not far from unitarity limit). In such a kinematic 

region a dipole of given transverse size b 2 expands with decrease of a; to a 

soft hadronic scale leading to a switch from the perturbative QCD unitarity 

regime to the soft QCD regime. A distinctive feature of such a scenario 

will be a fast increase with energy of the slope of the f-dependence of hard 

exclusive processes. 

A similar analysis of the black limit also applies to exclusive process such 

as diffractive photoproduction and electroproduction of J/ip. In this case, 

the black limit applies to the forward scattering amplitude: 

2itr 2 N / dzd 2 btp 1 (z,b 2 ,Q 2 )'il)j/^(z,b 



is 



(29) 



One can then further refine this, at a given x, by breaking the amplitude 
down into black, not-yet black and peripheral contributions in direct analogy 
with the decomposition for total cross sections in eq.(p6[). It follows from 
this formulae that within the black limit approximation the cross sections 
of exclusive processes at central impact parameters should not depend on 
energy. The difference from the structure functions case is because the 
wavefunction of vector meson suppresses the contribution of small transverse 
sizes. However, a residual, relatively slow, increase of amplitude with energy 
is expected due to contributions from peripheral collisions neglected in the 
above formulae. Since the analysis of quark Fermi motion effects Q ||] 
revealed that J / ip photoproduction is dominated by relatively large b, where 
the unitarity corrections set in early, this is an excellent process in which to 
search for the onset of the black limit. 
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Finally, proximity to Gribov's black body limit has important implica- 
tions for diffraction, which in this limit is equivalent to elastic scattering of 
quark-gluon configurations in the photon wave function and contributes half 
of the total cross section (recall that in the black limit a e i = Oi ne x = ctot/2). 
Since all component states in the photon have the same cross section non- 
diagonal transitions in M 2 are absent. We can therefore unfold the M 2 - 
integral in eq.(^7|) and divide by two to get to the diffractive mass spectrum: 

dF° {2 \x,Q 2 ) = nr 2 arget Q 2 M 2 p(M 2 ) 

dM 2 12tt 3 (M 2 + Q 2 ) 2 { ' 

Hence for M 2 ^> Q 2 one obtains an expression similar to the triple Pomeron 
limit with ajp(O) = 1. The contribution from scattering at peripheral im- 
pact parameters, neglected in the black body scenario, would correspond to 
Ojp(0) > 1. The distinctive feature of the final state is that the average 
transverse momenta are oc M/2. Let us assume that the black body limit is 
reached in nature, for example in DIS off heavy nuclei. Since it should be 
reached at lower energies for large dipoles, the diffractive production of small 
masses (M ~ few GeV) may serve as an early signal for the onset of the 
Gribov regime. Thus, for a given Q 2 eq.pO|) will be valid up to M 2 ~ Q 2 at 
so that decreasing x expands the region of applicability of this equation. 

Note also, as we have argued in that the current diffractive parton 
density analyses of the HERA diffractive data (see e.g [fl7f| ) indicate that 
gluon induced diffraction at HERA energies maybe already be close to this 
saturation limit. Unfortunately the lack of measurements of the t-slope 
of the diffractive amplitude around Q 2 ~ 4 — 6 GeV 2 in gluon induced 
channels precludes distinguishing between two competing scenarios. Either 
small sized colour octet dipoles with a e fj ~ 30-40 mb dominate (expected 
B ~ Bj/^ w 4.5 GeV -2 , a <C a so j t ,ajp — 1 ~ 0.2) or large colour-triplet 
dipole with with a e ff ~ 50-60 mb dominate (expected B ~ B inc i usive 7 
GeV -2 , ol ~ a so f t ,a]p - 1 ~ a so ft ). The first scenario naturally leads 
to a behaviour analogous to eq.(|30|) for the gluon channel for scattering off 
nuclei. We will consider these phenomena elsewhere. 

7 Conclusions 

We have proposed a physically motivated ansatz for the dipole cross section 
(DCS) relevant to a wide range of small x scattering processes. The small 
dipole cross section is governed by the leading log gluon density at small 
x. Using this and the measured pion-proton cross section as a guide, we 
construct an ansatz for it in the non-perturbative region, below the input 
scale at which the input density for the gluon is defined. At very small 
values of result of the large and steeply rising gluon density, the DCS 
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threatens to become larger at small perturbative b than the pion-proton 
cross section, in conflict with unitarity. To prevent this from happening we 
tame the rapid growth using a smooth ansatz that ensures a monotonically 
increasing function of b at fixed x. 

The resultant DCS produces values for Fl(x,Q 2 ) which are in good 
agreement with those from perturbative QCD in the large Q 2 and the high 
end of the small x region, where it would be expected to. Our DCS compares 
reasonably well with all available small x data on F2 from HERA, without 
any further tuning of parameters. Interestingly, in the moderate Q 2 region a 
significant fraction of the cross section appears to be coming from the region 
of large non-perturbative dipoles. If the perturbative unitarity corrections 
are neglected our model would continue to grow very steeply in the very 
small x region. More detailed studies of the choices that we make for the 
precise form of the ansatz are being carried out and will be reported in a 
separate paper [^(J. Although these clearly affect the results quantitatively 
to a certain extent, we are confident of the qualitative conclusion of the 
paper that unitarity corrections must set in within or near to the HERA 
kinematic region of small x and moderate Q 2 (we estimate that the Q 2 < 10 
GeV 2 will be affected at the smallest attainable values of x). This calls into 
question the use of the low, and x-independent, input scales used in the 
standard DGLAP fits % ^§ g|, H (Qo ~ °- 8 ~ 2 - 6 Gey2 ) 



In section(6) we generalized Gribov's black limit to the perturbative 
QCD analysis of DIS to make it compatible with the applicability of the 
DGLAP approximation at sufficiently large Q 2 but fixed x and with the 
dominance of peripheral collisions at fixed Q 2 at extremely high energies. 
We estimate that this black limit is close to the lower edge of the HERA 
kinematic region in x for Q 2 = 1 GeV 2 . We suggest that certain diffractive 
processes may acts as earlier indicators of the onset of this new regime. 
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